Multi-ancestry genome-wide meta-analysis identifies novel basal cell carcinoma loci and shared genetic effects with squamous cell carcinoma.

Basal cell carcinoma (BCC) is one of the most common malignancies worldwide, yet its genetic determinants are incompletely defined. We perform a European ancestry genome-wide association (GWA) meta-analysis and a Hispanic/Latino ancestry GWA meta-analysis and meta-analyze both in a multi-ancestry GWAS meta-analysis of BCC, totaling 50,531 BCC cases and 762,234 controls from four cohorts (GERA, Mass-General Brigham Biobank, UK Biobank, and 23andMe research cohort). Here we identify 122 BCC-associated loci, of which 36 were novel, and subsequently fine-mapped these associations. We also identify an association of the well-known pigment gene SLC45A2 as well as associations at RCC2 and CLPTM1L with BCC in Hispanic/Latinos. We examine these BCC loci for association with cutaneous squamous cell carcinoma (cSCC) in 16,407 SCC cases and 762,486 controls of European ancestry, and 33 SNPs show evidence of association. Our study findings provide important insights into the genetic basis of BCC and cSCC susceptibility.

European and multi-ancestry genome-wide association meta-analysis of atopic dermatitis highlights importance of systemic immune regulation.

Atopic dermatitis (AD) is a common inflammatory skin condition and prior genome-wide association studies (GWAS) have identified 71 associated loci. In the current study we conducted the largest AD GWAS to date (discovery N = 1,086,394, replication N = 3,604,027), combining previously reported cohorts with additional available data. We identified 81 loci (29 novel) in the European-only analysis (which all replicated in a separate European analysis) and 10 additional loci in the multi-ancestry analysis (3 novel). Eight variants from the multi-ancestry analysis replicated in at least one of the populations tested (European, Latino or African), while two may be specific to individuals of Japanese ancestry. AD loci showed enrichment for DNAse I hypersensitivity and eQTL associations in blood. At each locus we prioritised candidate genes by integrating multi-omic data. The implicated genes are predominantly in immune pathways of relevance to atopic inflammation and some offer drug repurposing opportunities.

Genotyping, sequencing and analysis of 140,000 adults from Mexico City.

The Mexico City Prospective Study is a prospective cohort of more than 150,000 adults recruited two decades ago from the urban districts of Coyoacán and Iztapalapa in Mexico City1. Here we generated genotype and exome-sequencing data for all individuals and whole-genome sequencing data for 9,950 selected individuals. We describe high levels of relatedness and substantial heterogeneity in ancestry composition across individuals. Most sequenced individuals had admixed Indigenous American, European and African ancestry, with extensive admixture from Indigenous populations in central, southern and southeastern Mexico. Indigenous Mexican segments of the genome had lower levels of coding variation but an excess of homozygous loss-of-function variants compared with segments of African and European origin. We estimated ancestry-specific allele frequencies at 142 million genomic variants, with an effective sample size of 91,856 for Indigenous Mexican ancestry at exome variants, all available through a public browser. Using whole-genome sequencing, we developed an imputation reference panel that outperforms existing panels at common variants in individuals with high proportions of central, southern and southeastern Indigenous Mexican ancestry. Our work illustrates the value of genetic studies in diverse populations and provides foundational imputation and allele frequency resources for future genetic studies in Mexico and in the United States, where the Hispanic/Latino population is predominantly of Mexican descent.

Association of Behavioral and Clinical Risk Factors With Cataract: A Two-Sample Mendelian Randomization Study.

Purpose: To investigate the association of genetically determined primary open-angle glaucoma (POAG), myopic refractive error (RE), type 2 diabetes (T2D), blood pressure (BP), body mass index (BMI), cigarette smoking, and alcohol consumption with the risk of age-related cataract. Methods: To assess potential causal effects of clinical or behavioral factors on cataract risk, we conducted two-sample Mendelian randomization analyses. Genetic instruments, based on common genetic variants associated with risk factors at genome-wide significance (P < 5 × 10-8), were derived from published genome-wide association studies (GWAS). For age-related cataract, we used GWAS summary statistics from our previous GWAS conducted in the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort (28,092 cataract cases and 50,487 controls; all non-Hispanic whites) or in the UK Biobank (31,852 cataract cases and 428,084 controls; all European-descent individuals). We used the inverse-variance weighted (IVW) method as our primary source of Mendelian randomization estimates and conducted common sensitivity analyses. Results: We found that genetically determined POAG and mean spherical equivalent RE were significantly associated with cataract risk (IVW model: odds ratio [OR] = 1.04; 95% confidence interval [CI], 1.01-1.08; P = 0.018; per diopter more hyperopic: OR = 0.92; 95% CI, 0.89-0.93; P = 6.51 × 10-13, respectively). In contrast, genetically determined T2D, BP, BMI, cigarette smoking, or alcohol consumption were not associated with cataract risk (P > 0.05). Conclusions: Our results provide evidence that genetic risks for POAG and myopia may be causal risk factors for age-related cataract. These results are consistent with previous observational studies reporting associations of myopia with cataract risk. This information may support population cataract risk stratification and screening strategies.

Rare coding variants in CHRNB2 reduce the likelihood of smoking.

Human genetic studies of smoking behavior have been thus far largely limited to common variants. Studying rare coding variants has the potential to identify drug targets. We performed an exome-wide association study of smoking phenotypes in up to 749,459 individuals and discovered a protective association in CHRNB2, encoding the ß2 subunit of the α4ß2 nicotine acetylcholine receptor. Rare predicted loss-of-function and likely deleterious missense variants in CHRNB2 in aggregate were associated with a 35% decreased odds for smoking heavily (odds ratio (OR) = 0.65, confidence interval (CI) = 0.56-0.76, P = 1.9 × 10-8). An independent common variant association in the protective direction ( rs2072659 ; OR = 0.96; CI = 0.94-0.98; P = 5.3 × 10-6) was also evident, suggesting an allelic series. Our findings in humans align with decades-old experimental observations in mice that ß2 loss abolishes nicotine-mediated neuronal responses and attenuates nicotine self-administration. Our genetic discovery will inspire future drug designs targeting CHRNB2 in the brain for the treatment of nicotine addiction.

Large-scale multitrait genome-wide association analyses identify hundreds of glaucoma risk loci.

Glaucoma, a leading cause of irreversible blindness, is a highly heritable human disease. Previous genome-wide association studies have identified over 100 loci for the most common form, primary open-angle glaucoma. Two key glaucoma-associated traits also show high heritability: intraocular pressure and optic nerve head excavation damage quantified as the vertical cup-to-disc ratio. Here, since much of glaucoma heritability remains unexplained, we conducted a large-scale multitrait genome-wide association study in participants of European ancestry combining primary open-angle glaucoma and its two associated traits (total sample size over 600,000) to substantially improve genetic discovery power (263 loci). We further increased our power by then employing a multiancestry approach, which increased the number of independent risk loci to 312, with the vast majority replicating in a large independent cohort from 23andMe, Inc. (total sample size over 2.8 million; 296 loci replicated at P < 0.05, 240 after Bonferroni correction). Leveraging multiomics datasets, we identified many potential druggable genes, including neuro-protection targets likely to act via the optic nerve, a key advance for glaucoma because all existing drugs only target intraocular pressure. We further used Mendelian randomization and genetic correlation-based approaches to identify novel links to other complex traits, including immune-related diseases such as multiple sclerosis and systemic lupus erythematosus.

A multiethnic genome-wide analysis of 19,420 individuals identifies novel loci associated with axial length and shared genetic influences with refractive error and myopia.

Introduction: Long axial length (AL) is a risk factor for myopia. Although family studies indicate that AL has an important genetic component with heritability estimates up to 0.94, there have been few reports of AL-associated loci. Methods: Here, we conducted a multiethnic genome-wide association study (GWAS) of AL in 19,420 adults of European, Latino, Asian, and African ancestry from the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort, with replication in a subset of the Consortium for Refractive Error and Myopia (CREAM) cohorts of European or Asian ancestry. We further examined the effect of the identified loci on the mean spherical equivalent (MSE) within the GERA cohort. We also performed genome-wide genetic correlation analyses to quantify the genetic overlap between AL and MSE or myopia risk in the GERA European ancestry sample. Results: Our multiethnic GWA analysis of AL identified a total of 16 genomic loci, of which 5 are novel. We found that all AL-associated loci were significantly associated with MSE after Bonferroni correction. We also found that AL was genetically correlated with MSE (rg = -0.83; SE, 0.04; p = 1.95 × 10-89) and myopia (rg = 0.80; SE, 0.05; p = 2.84 × 10-55). Finally, we estimated the array heritability for AL in the GERA European ancestry sample using LD score regression, and found an overall heritability estimate of 0.37 (s.e. = 0.04). Discussion: In this large and multiethnic study, we identified novel loci, associated with AL at a genome-wide significance level, increasing substantially our understanding of the etiology of AL variation. Our results also demonstrate an association between AL-associated loci and MSE and a shared genetic basis between AL and myopia risk.

MiXcan: a framework for cell-type-aware transcriptome-wide association studies with an application to breast cancer.

Human bulk tissue samples comprise multiple cell types with diverse roles in disease etiology. Conventional transcriptome-wide association study approaches predict genetically regulated gene expression at the tissue level, without considering cell-type heterogeneity, and test associations of predicted tissue-level expression with disease. Here we develop MiXcan, a cell-type-aware transcriptome-wide association study approach that predicts cell-type-level expression, identifies disease-associated genes via combination of cell-type-level association signals for multiple cell types, and provides insight into the disease-critical cell type. As a proof of concept, we conducted cell-type-aware analyses of breast cancer in 58,648 women and identified 12 transcriptome-wide significant genes using MiXcan compared with only eight genes using conventional approaches. Importantly, MiXcan identified genes with distinct associations in mammary epithelial versus stromal cells, including three new breast cancer susceptibility genes. These findings demonstrate that cell-type-aware transcriptome-wide analyses can reveal new insights into the genetic and cellular etiology of breast cancer and other diseases.

Common and rare variant associations with clonal haematopoiesis phenotypes.

Clonal haematopoiesis involves the expansion of certain blood cell lineages and has been associated with ageing and adverse health outcomes1-5. Here we use exome sequence data on 628,388 individuals to identify 40,208 carriers of clonal haematopoiesis of indeterminate potential (CHIP). Using genome-wide and exome-wide association analyses, we identify 24 loci (21 of which are novel) where germline genetic variation influences predisposition to CHIP, including missense variants in the lymphocytic antigen coding gene LY75, which are associated with reduced incidence of CHIP. We also identify novel rare variant associations with clonal haematopoiesis and telomere length. Analysis of 5,041 health traits from the UK Biobank (UKB) found relationships between CHIP and severe COVID-19 outcomes, cardiovascular disease, haematologic traits, malignancy, smoking, obesity, infection and all-cause mortality. Longitudinal and Mendelian randomization analyses revealed that CHIP is associated with solid cancers, including non-melanoma skin cancer and lung cancer, and that CHIP linked to DNMT3A is associated with the subsequent development of myeloid but not lymphoid leukaemias. Additionally, contrary to previous findings from the initial 50,000 UKB exomes6, our results in the full sample do not support a role for IL-6 inhibition in reducing the risk of cardiovascular disease among CHIP carriers. Our findings demonstrate that CHIP represents a complex set of heterogeneous phenotypes with shared and unique germline genetic causes and varied clinical implications.

Assessment of genetic susceptibility to multiple primary cancers through whole-exome sequencing in two large multi-ancestry studies.

BACKGROUND: Up to one of every six individuals diagnosed with one cancer will be diagnosed with a second primary cancer in their lifetime. Genetic factors contributing to the development of multiple primary cancers, beyond known cancer syndromes, have been underexplored. METHODS: To characterize genetic susceptibility to multiple cancers, we conducted a pan-cancer, whole-exome sequencing study of individuals drawn from two large multi-ancestry populations (6429 cases, 165,853 controls). We created two groupings of individuals diagnosed with multiple primary cancers: (1) an overall combined set with at least two cancers across any of 36 organ sites and (2) cancer-specific sets defined by an index cancer at one of 16 organ sites with at least 50 cases from each study population. We then investigated whether variants identified from exome sequencing were associated with these sets of multiple cancer cases in comparison to individuals with one and, separately, no cancers. RESULTS: We identified 22 variant-phenotype associations, 10 of which have not been previously discovered and were significantly overrepresented among individuals with multiple cancers, compared to those with a single cancer. CONCLUSIONS: Overall, we describe variants and genes that may play a fundamental role in the development of multiple primary cancers and improve our understanding of shared mechanisms underlying carcinogenesis.

ANGPTL7, a therapeutic target for increased intraocular pressure and glaucoma.

Glaucoma is a leading cause of blindness. Current glaucoma medications work by lowering intraocular pressure (IOP), a risk factor for glaucoma, but most treatments do not directly target the pathological changes leading to increased IOP, which can manifest as medication resistance as disease progresses. To identify physiological modulators of IOP, we performed genome- and exome-wide association analysis in >129,000 individuals with IOP measurements and extended these findings to an analysis of glaucoma risk. We report the identification and functional characterization of rare coding variants (including loss-of-function variants) in ANGPTL7 associated with reduction in IOP and glaucoma protection. We validated the human genetics findings in mice by establishing that Angptl7 knockout mice have lower (~2 mmHg) basal IOP compared to wild-type, with a trend towards lower IOP also in heterozygotes. Conversely, increasing murine Angptl7 levels via injection into mouse eyes increases the IOP. We also show that acute Angptl7 silencing in adult mice lowers the IOP (~2-4 mmHg), reproducing the observations in knockout mice. Collectively, our data suggest that ANGPTL7 is important for IOP homeostasis and is amenable to therapeutic modulation to help maintain a healthy IOP that can prevent onset or slow the progression of glaucoma.

Association Between Myopic Refractive Error and Primary Open-Angle Glaucoma: A 2-Sample Mendelian Randomization Study.

Importance: Refractive error (RE) is the most common form of visual impairment, and myopic RE is associated with an increased risk of primary open-angle glaucoma (POAG). Whether this association represents a causal role of RE in the etiology of POAG remains unknown. Objective: To evaluate shared genetic influences and investigate the association of myopic RE with the risk for POAG. Design, Setting, and Participants: Observational analyses were used to evaluate the association between mean spherical equivalent (MSE) RE (continuous trait) or myopia (binary trait) and POAG risk in individuals from the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort. To quantify genetic overlap, genome-wide genetic correlation analyses were performed using genome-wide association studies (GWAS) of MSE RE or myopia and POAG from GERA. Potential causal effects were assessed between MSE RE and POAG using 2-sample Mendelian randomization. Genetic variants associated with MSE RE were derived using GWAS summary statistics from a GWAS of RE conducted in 102 117 UK Biobank participants. For POAG, we used GWAS summary statistics from our previous GWAS (3836 POAG cases and 48 065 controls from GERA). Data analyses occurred between July 2020 and October 2021. Main Outcomes and Measure: Our main outcome was POAG risk as odds ratio (OR) caused by per-unit difference in MSE RE (in diopters). Results: Our observational analyses included data for 54 755 non-Hispanic White individuals (31 926 [58%] females and 22 829 [42%] males). Among 4047 individuals with POAG, mean (SD) age was 73.64 (9.20) years; mean (SD) age of the 50 708 controls was 65.38 (12.24) years. Individuals with POAG had a lower refractive MSE and were more likely to have myopia or high myopia compared with the control participants (40.2% vs 34.1%, P = 1.31 × 10-11 for myopia; 8.5% vs 6.8%, P = .004 for high myopia). Our genetic correlation analyses demonstrated that POAG was genetically correlated with MSE RE (rg, -0.24; SE, 0.06; P = 3.90 × 10-5), myopia (rg, 0.21; SE, 0.07; P = .004), and high myopia (rg, 0.23; SE, 0.09; P = .01). Genetically assessed refractive MSE was negatively associated with POAG risk (inverse-variance weighted model: OR per diopter more hyperopic MSE = 0.94; 95% CI, 0.89-0.99; P = .01). Conclusions and Relevance: These findings demonstrate a shared genetic basis and an association between myopic RE and POAG risk. This may support population POAG risk stratification and screening strategies, based on RE information.

Genome-wide association study of actinic keratosis identifies new susceptibility loci implicated in pigmentation and immune regulation pathways.

Actinic keratosis (AK) is a common precancerous cutaneous neoplasm that arises on chronically sun-exposed skin. AK susceptibility has a moderate genetic component, and although a few susceptibility loci have been identified, including IRF4, TYR, and MC1R, additional loci have yet to be discovered. We conducted a genome-wide association study of AK in non-Hispanic white participants of the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort (n = 63,110, discovery cohort), with validation in the Mass-General Brigham (MGB) Biobank cohort (n = 29,130). We identified eleven loci (P < 5 × 10-8), including seven novel loci, of which four novel loci were validated. In a meta-analysis (GERA + MGB), one additional novel locus, TRPS1, was identified. Genes within the identified loci are implicated in pigmentation (SLC45A2, IRF4, BNC2, TYR, DEF8, RALY, HERC2, and TRPS1), immune regulation (FOXP1 and HLA-DQA1), and cell signaling and tissue remodeling (MMP24) pathways. Our findings provide novel insight into the genetics and pathogenesis of AK susceptibility.

Genome-Wide Association Study Identifies Two Common Loci Associated with Pigment Dispersion Syndrome/Pigmentary Glaucoma and Implicates Myopia in its Development.

PURPOSE: To identify genetic variants associated with pigment dispersion syndrome (PDS) and pigmentary glaucoma (PG) in unrelated patients and to further understand the genetic and potentially causal relationships between PDS and associated risk factors. DESIGN: A 2-stage genome-wide association meta-analysis with replication and subsequent in silico analyses including Mendelian randomization. PARTICIPANTS: A total of 574 cases with PG or PDS and 52 627 controls of European descent. METHODS: Genome-wide association analyses were performed in 4 cohorts and meta-analyzed in 3 stages: (1) a discovery meta-analysis was performed in 3 cohorts, (2) replication was performed in the fourth cohort, and (3) all 4 cohorts were meta-analyzed to increase statistical power. Two-sample Mendelian randomization was used to determine whether refractive error and intraocular pressure exert causal effects over PDS. MAIN OUTCOME MEASURES: The association of genetic variants with PDS and whether myopia exerts causal effects over PDS. RESULTS: Significant association was present at 2 novel loci for PDS/PG. These loci and follow-up analyses implicate the genes gamma secretase activator protein (GSAP) (lead single nucleotide polymorphism [SNP]: rs9641220, P = 6.0×10-10) and glutamate metabotropic receptor 5 (GRM5)/TYR (lead SNP: rs661177, P = 3.9×10-9) as important factors in disease risk. Mendelian randomization showed significant evidence that negative refractive error (myopia) exerts a direct causal effect over PDS (P = 8.86×10-7). CONCLUSIONS: Common SNPs relating to the GSAP and GRM5/TYR genes are associated risk factors for the development of PDS and PG. Although myopia is a known risk factor, this study uses genetic data to demonstrate that myopia is, in part, a cause of PDS and PG.

Ancestry- and sex-specific effects underlying inguinal hernia susceptibility identified in a multiethnic genome-wide association study meta-analysis.

Inguinal hernias are some of the most frequently diagnosed conditions in clinical practice and inguinal hernia repair is the most common procedure performed by general surgeons. Studies of inguinal hernias in non-European populations are lacking, though it is expected that such studies could identify novel loci. Further, the cumulative lifetime incidence of inguinal hernia is nine times greater in men than women, however, it is not clear why this difference exists. We conducted a genome-wide association meta-analysis of inguinal hernia risk across 513 120 individuals (35 774 cases and 477 346 controls) of Hispanic/Latino, African, Asian and European descent, with replication in 728 418 participants (33 491 cases and 694 927 controls) from the 23andMe, Inc dataset. We identified 63 genome-wide significant loci (P < 5 × 10-8), including 41 novel. Ancestry-specific analyses identified two loci (LYPLAL1-AS1/SLC30A10 and STXBP6-NOVA1) in African ancestry individuals. Sex-stratified analyses identified two loci (MYO1D and ZBTB7C) that are specific to women, and four (EBF2, EMX2/RAB11FIP2, VCL and FAM9A/FAM9B) that are specific to men. Functional experiments demonstrated that several of the associated regions (EFEMP1 and LYPLAL1-SLC30A10) function as enhancers and show differential activity between risk and reference alleles. Our study highlights the importance of large-scale genomic studies in ancestrally diverse populations for identifying ancestry-specific inguinal hernia susceptibility loci and provides novel biological insights into inguinal hernia etiology.

Whole genome sequence analysis of platelet traits in the NHLBI Trans-Omics for Precision Medicine (TOPMed) initiative.

Platelets play a key role in thrombosis and hemostasis. Platelet count (PLT) and mean platelet volume (MPV) are highly heritable quantitative traits, with hundreds of genetic signals previously identified, mostly in European ancestry populations. We here utilize whole genome sequencing (WGS) from NHLBI's Trans-Omics for Precision Medicine initiative (TOPMed) in a large multi-ethnic sample to further explore common and rare variation contributing to PLT (n = 61 200) and MPV (n = 23 485). We identified and replicated secondary signals at MPL (rs532784633) and PECAM1 (rs73345162), both more common in African ancestry populations. We also observed rare variation in Mendelian platelet-related disorder genes influencing variation in platelet traits in TOPMed cohorts (not enriched for blood disorders). For example, association of GP9 with lower PLT and higher MPV was partly driven by a pathogenic Bernard-Soulier syndrome variant (rs5030764, p.Asn61Ser), and the signals at TUBB1 and CD36 were partly driven by loss of function variants not annotated as pathogenic in ClinVar (rs199948010 and rs571975065). However, residual signal remained for these gene-based signals after adjusting for lead variants, suggesting that additional variants in Mendelian genes with impacts in general population cohorts remain to be identified. Gene-based signals were also identified at several genome-wide association study identified loci for genes not annotated for Mendelian platelet disorders (PTPRH, TET2, CHEK2), with somatic variation driving the result at TET2. These results highlight the value of WGS in populations of diverse genetic ancestry to identify novel regulatory and coding signals, even for well-studied traits like platelet traits.

A large-scale transcriptome-wide association study (TWAS) of 10 blood cell phenotypes reveals complexities of TWAS fine-mapping.

Hematological measures are important intermediate clinical phenotypes for many acute and chronic diseases and are highly heritable. Although genome-wide association studies (GWAS) have identified thousands of loci containing trait-associated variants, the causal genes underlying these associations are often uncertain. To better understand the underlying genetic regulatory mechanisms, we performed a transcriptome-wide association study (TWAS) to systematically investigate the association between genetically predicted gene expression and hematological measures in 54,542 Europeans from the Genetic Epidemiology Research on Aging cohort. We found 239 significant gene-trait associations with hematological measures; we replicated 71 associations at p < 0.05 in a TWAS meta-analysis consisting of up to 35,900 Europeans from the Women's Health Initiative, Atherosclerosis Risk in Communities Study, and BioMe Biobank. Additionally, we attempted to refine this list of candidate genes by performing conditional analyses, adjusting for individual variants previously associated with hematological measures, and performed further fine-mapping of TWAS loci. To facilitate interpretation of our findings, we designed an R Shiny application to interactively visualize our TWAS results by integrating them with additional genetic data sources (GWAS, TWAS from multiple reference panels, conditional analyses, known GWAS variants, etc.). Our results and application highlight frequently overlooked TWAS challenges and illustrate the complexity of TWAS fine-mapping.